The r-mode instability in strange stars with a crystalline crust 
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The r-mode instability, believed to limit the rotation speed of compact stars, can provide empirical 
confirmation for the existence of stable deconfined phases of quark matter that are predicted by weak 
coupling calculations in Quantum Chromodynamics. We construct a model for strange quark stars 
as heavy as 2Mq that are made of superconducting quark matter in the bulk and a thin crystalline 
quark matter crust. This crystalline quark crust is sufficiently robust to withstand r-mode heating 
and viscous rubbing for realistic mode amplitudes 0(1O -2 ), unlike a crust made of neutron-rich 
nuclei. The dissipation provided by viscous rubbing at the core-crust boundary is both necessary 
and sufficient to obtain stable rotation speeds that are consistent with the majority of rapidly 
spinning pulsars in low mass X-ray binaries. Our analysis implies that while bare strange stars are 
ruled out by the existence of rapidly spinning pulsars, a strange star with a quark matter crust is a 
distinct possibility. 
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Introduction. — Compact stars provide a unique win- 
dow into the physics of matter at densities well above nu- 
clear saturation. At asymptotically high densities, Quan- 
tum Chromodynamics (QCD) favors a maximally sym- 
metric phase of superconducting quark matter called the 
color-flavor-locked (CFL) phase [T]. However, at com- 
pact star densities, the phase of matter is less certain. In 
this intermediate density regime, some of the most im- 
portant constraints on theoretical models of dense matter 
come from compact star observables, eg., mass-radius re- 
lations [2] and cooling rates in compact stars [3]. The 
spin frequencies of fast spinning pulsars is yet another, 
but herein lies a puzzle. 

The spin frequency v 716 Hz of the fastest spin- 
ning neutron star PSR J1748-2446ad [I] is still well be- 
low the estimated theoretical maximum « 1000 Hz cor- 
responding to the break-up Kepler frequency. The ma- 
jority of "spun-up" neutron stars in low mass X-ray bi- 
naries (LMXBs) cluster between 300-640 Hz [5]. What 
determines the maximum spin frequency and how does it 
relate to the nature of dense matter inside the star? It is 
believed that the r-mode instability [HE] sets the limit on 
spin frequency. The r-mode is driven unstable by gravita- 
tional wave emission at the expense of rotational energy, 
thereby spinning the star down. Internal sources of fric- 
tion such as bulk/shear viscosity tend to damp r-mode 
growth. At any given temperature T, the competition be- 
tween the former growing on a characteristic time scale 
t gw and the latter acting on a time scale rp determines 
the limiting or critical frequency f2 c . The gravitational 
waves from r-mode oscillations may be detected given the 
expected improved sensitivities of ground-based interfer- 
ometers (eg. VIRGO, Advanced LIGO 012]), providing 
a novel diagnostic of dense matter in compact stars. 



In this letter, we investigate the maximum spin fre- 
quency of strange stars that can sustain the r-mode in- 
stability. If discovered, strange stars, which are believed 
to be made entirely of deconfined quarks, would confirm 
Witten's hypothesis of self-bound strange quark mat- 
ter [TO] as the true ground state of matter. At the theo- 
retical level, the r-mode mechanism for limiting the star's 
spin frequency does not seem to work for bare strange 
stars, i.e, one without a crust pTHTB"] . The critical fre- 
quency falls short of typical LMXB frequencies, essen- 
tially because of insufficient dissipation of the r-mode en- 
ergy. However, studies of ordinary neutron stars |16H18j 
suggest that r-modes can be strongly affected (damped) 
by the presence of a nuclear crust. This extra dissipation 
compared to a pure fluid is mediated by kinetic viscos- 
ity in a thin boundary layer and significantly shortens 
rp. Consequently, the critical frequency £l c increases dra- 
matically. Such a large effect is not expected for neutron 
stars because even r-modes with relatively small ampli- 
tude O(10~ 2 ) can melt the nuclear crust [17J . But com- 
pared to a nuclear crust, the crystalline quark phase is 
much more rigid, with a shear modulus that is 20 — 1000 
larger [191 and a melting temperature Tmdt ~ A w 10 
MeV that is 10 times larger. If crystalline quark matter, 
rather than nuclei, is taken to constitute the crust, it al- 
lows for a large damping of the r-mode that brings the 
strange star model into agreement with observed LMXB 
spin rates. We find that if one is to entertain the Witten 
hypothesis of the absolute stability of quark matter at 
high density, self-bound strange stars likely have a crust 
made of superconducting quarks in the crystalline phase. 

Self-bound stars with a crystalline crust. — At realis- 
tic compact star densities with quark chemical potential 
jj, q ~ 300 MeV, the strange quark mass m s ~ 150 MeV is 
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non-negligible, and severely stresses the CFL phase. The 
energetically favored phase is a chirally rotated phase of 
CFL [20] which leads to kaon condensation when the en- 
ergy cost m 2 8 /{2p q ) exceeds the CFL kaon mass tuk- A 
homogeneous strange star in this kaon condensed CFL 
phase (henceforth CFL-KO) has already been ruled out 
from an r-mode study [TS]. We consider the effect of a 
thin layer of superconducting quark matter in the crys- 
talline phase forming a crust on the surface of a CFL-KO 
star. In the CFL-KO phase, the equation of state can be 
written as [2U 
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where B is the MIT Bag constant, the decay constant 
fl = (21 - 81n2)/(367r 2 ), and a 4 ~ 0.7 is a reasonable 
choice for perturbative QCD-inspired corrections [22] . 
This parameterization allows for quark stars as heavy as 
2Mq and a phase transition to quark matter above nu- 
clear saturation. For a fixed value of the strange quark 
mass m s , CFL-KO gap parameter A and 04, the bag con- 
stant B is determined by our choice of a given binding en- 



ergy per baryon E/A for bulk strange quark matter. The 
mass M and radius R of the compact star is calculated us- 
ing the Tolman-Oppenheimer-Volkov equations with gen- 
eral relativistic corrections. For fixed m s , A, 04 and B, 
the central density p c of the star is adjusted as a factor of 
the nuclear saturation density po = 2.55 x 10 14 g/cm 3 to 
produce different stellar masses. The core-crust bound- 
ary is determined by comparing the free energy of the 
CFL-KO with the most favored crystalline quark phase at 
these densities 23 . This is the CubeX phase, a specific 
generalization of the LOFF phase where the quark pairs 
carry non-zero total momentum 2q n [TU] . The gap varies 
in space as ~ exp(iq„ • x) according to the reciprocal 
vectors {q n } which define the specific crystal structure of 
this phase (as opposed to a physically rigid lattice) [2~4"] . 
The chemical potential at which the crystalline phase 
takes over is quite sensitive to the gap, hence we tune 
the gap to locate this transition point near the surface of 
the star at radius r c < R. Though the gap, and hence 
the location of the crust cannot be known for certain, 
the above construction provides the most instructive sce- 
nario, in that it allows us to demonstrate the effect that 
even a thin crust d crU st = R — r c <§; R can have in provid- 
ing additional damping compared to homogeneous quark 
phases. Table [I] shows the two sets of parameters we use 
in our calculations. 



m a (MeV) 
160 
160 



B 1/4 (MeV) 
137.8 
136.1 



E/A (MeV) 
930 
920 



Q.4 

0.7 
0.7 



A (MeV) 
24.0 
24.2 



Pc(po) 
2.71 
4.93 



R (km) 
11.14 
11.71 



M(M ) 
1.40 
1.97 



MR) 
0.997 
0.993 



fi{r c ) (MeV) 
310.3 
307.7 



H(R) (MeV) 
310.0 
306.7 



TABLE I. Parameters for two self-bound stars, M/M©=1.40 and 1.97. The location of the crust in stellar units is r c , whereat 
the transition chemical potential is fi(r c ). 



r-mode damping and critical spin frequency. — The 
restoring force for the r-mode oscillation is the Coriolis 
force, and ignoring viscous effects, the mode frequency 
uj r in the incrtial frame is [7] 
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where uj mt is the mode frequency in the frame that co- 
rotates with the star, and m=2 is the azimuthal quantum 
number of the first r-mode to become unstable. Higher 
order corrections include the sensitivity to the density 
profile p(r) (see Refs. [13, 25]). The r-mode energy grows 
at the expense of the star's rotational kinetic energy as 
angular momentum is lost to gravitational wave emission. 
In general, viscous forces counter this energy growth. 
The time scale r associated with growth or dissipation 



is given by 1/Tj = — ^ (^f)-> where the subscript cor- 
responds to the growth or dissipative mechanisms and 
E is the mode energy. The time scale for gravitational 
radiation is given by [7] 
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This timcscale is negative, indicating exponential mode 
growth (the instability). The damping timescales from 
shear viscosity t v and bulk viscosity for the CFL-KO 
phase have been calculated in Ref. [15]. The important 
new contribution comes from the damping in the vis- 
cous boundary layer (VBL) just beneath the rigid crust 
located at r c (Table The r-mode oscillations decay 



in the crust exponentially within a characteristic length 
scale [17] 
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where rj{r c ) is the shear viscosity of the matter at den- 
sity p(r c ). For the range of parameters we explored, 
the crust thickness B — r c is always at least a factor of 
about 100 larger than d, so the VBL approximation is 
self-consistent. The damping time associated with dissi- 
pation in the VBL is called the rubbing time scale, given 

by m 
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where I m f=s 0.8041. The above estimate assumes that 
mechanical coupling between the core and the crust is 
absent, but if the r-mode oscillations penetrate into the 
crust, it increases the rubbing time scale by a "slippage" 
factor of (Su/u)~ 2 [THl [T7], where u is the oscillating 
fluid velocity in the core and 5u the difference in the 
fluid velocity in the core and the crust across the bound- 
ary. For a nuclear crust, the slippage factor spans a 
range 0.003 < (5u/u) 2 < 1 over a wide range of spin 
frequencies < Q/CIk < 0.5. Qk = %\T^Gj> is the 
Kepler frequency with p being the average matter den- 
sity. Crystalline quark matter is much more rigid than 
the nuclear crust. Performing a mode penetration cal- 
culation for the quark crust along the lines of Ref . [TB] , 
we find that the slippage factor for a canonical r-mode 
frequency lu — 20/3 is nearly 1 for the entire frequency 
range of interest, assuming a shear modulus about 100 
times larger than the nuclear crust. This is reasonable 
given how rigid the quark crust is. 

The critical frequency of the star at any temperature 
can be determined by the criterion that at this frequency, 
the fraction of energy dissipated per unit time by viscous 
processes exactly cancels the fraction of energy fed per 
unit time into the r-mode by gravitational wave emission: 
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The critical frequency curves obtained from solving 
Eq. ^ are plotted in Fig. [T] where at frequencies be- 
low the curve the r-mode oscillations are suppressed 
due to the dissipative processes. We choose param- 
eters for a small kaon condensation such that 5m = 
rriK — {m 2 — m 2 d )/{2p q ) = —2 MeV, with down quark 
mass rrid — 7 MeV. The kaon interaction coupling C is 
varied over a reasonable range 0.1 to 10 [2B]. Weaker 
coupling gives larger shear viscosity. The shear viscosity 
is sensitive to the kaon condensation, increasing sharply 



with decreasing 5m j^Hj- However, one cannot decrease 
the coupling C or 5m arbitrarily to increase the critical 
frequency and explain the LMXB spin rates. Too weak 
a coupling or small 5m would make the mean free path 
of the kaons larger than the star radius invalidating the 
hydrodynamic calculation. For example, the dip in the 
critical frequency around 7 x 10 9 K corresponds to ex- 
cluding phonon contributions when its mean free path 
becomes comparable to the star radius. We have shown 
that in the absence of the crust, even a small C is not 
sufficient to reproduce LMXB spin frequencies [TS]. We 
also limit the critical frequency to be below the limit set 
by the Kepler frequency VLk ■ This results in the plateaus 

near T = 10 10 K. We find in Fig. [l] that at 5m 2 

MeV, a small C ~ 0.1 is compatible with the data. 
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FIG. 1. Critical frequency for the M = M/M© = 1.97 star 
(blue curve) and M — M/Mq = 1.4 star (red curve). 

Including only the dissipation time from the VBL r ru b 
and equating it to the gravitational time scale tqw gives 
the critical frequency 
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where we used parameters for the 1.4M Q star with C — 
0.1. At T w 10 8 K, u c fa 350 Hz, thus we see that rubbing 
friction alone can generate sufficient damping to explain 
typical LMXB frequencies. In contrast, equating tqw 
to t v at T — 10 s K gives v c w 11 Hz. These estimates 
clearly highlight the essential role of the quark crust in 
providing a large damping of the r-mode. 

Discussion. — r-mode damping in a self-bound strange 
quark star with a crystalline quark crust is examined. 
For a plausible range of crust shear modulus, the canon- 
ical r-mode is damped at the core-crust boundary in 
exactly the right amount to tame the instability at ob- 
served LMXB frequencies and temperatures. Additional 
damping agents, such as mutual friction associated with 
the phonons and the kaon-condensed phase could pro- 
vide some extra small dissipation, but the critical spin 
frequency estimates presented here are a robust lower 
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bound. Further, we use an equation of state that meets 
the observational constraint of a large maximum mass of 
2Mq, making it a reasonable starting point for further 
investigation into properties of self-bound stars. Such 
stars most likely should have a thin crust made of super- 
conducting quarks in the crystalline phase if the Witten 
hypothesis of the absolute stability of quark matter at 
high density holds true. Finally, physical applications 
of our model might come from analyzing shear modes 
of the crystalline quark crust and matching them to fre- 
quencies extracted from seismic vibrations in magnetar 
flares, which appear to pose a challenge for normal quark 
matter crusts [27] . 
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